Identification of biomarker genes from multiple studies for abiotic stress in maize through machine learning.

Abiotic stresses are major limiting factors for maize growth. Therefore, exploration of the mechanisms underlying the response to abiotic stress in maize is of great interest. Toward this end, we performed integration of the feature selection method into the meta-analysis of microarray gene expression. Following extraction of raw data, normalization, and batch effect removal, the data were merged into one expression profile. Differentially expressed genes (DEGs) between control and abiotic conditions were used for the feature selection algorithm to find the minimum features for high-performance classification. Feature selection was performed using a correlation-based feature selection (CFS) algorithm, considering features with a coefficient of 0.7 to 1. Different algorithms of Bayes, Functions, Lazy, Meta, Rules, and Trees were then tested in order to classify the samples and find the best performance classifier in each group. Moreover, the biological pathways and promoter motif analysis of selected genes were identified. The superior and overall performance of classification using all features (DEGs) were 98.86% (Multilayer Perceptron) and 81.25%, respectively. Classification based on feature selection resulted in an average accuracy of 94.69% and 93.56% with 33 and 12 features, respectively. Subsequently, gene ontology and promoter analysis were performed for the 12 selected biomarker genes. Five of them were downregulated and 7 were upregulated. ABRE, unnamed-1, G-box, and G-Box are motifs related to genes involved in several abiotic stress responses and are located upstream of at least nine probes in our study. This study revealed key genes associated with tolerance to abiotic stress in maize.

A Genomic BSAseq Approach for the Characterization of QTLs Underlying Resistance to Fusarium oxysporum in Eggplant.

Eggplant (Solanum melongena L.), similar to many other crops, suffers from soil-borne diseases, including Fusarium oxysporum f. sp. melongenae (Fom), causing wilting and heavy yield loss. To date, the genetic factors underlying plant responses to Fom are not well known. We previously developed a Recombinant Inbred Lines (RILs) population using as a female parent the fully resistant line '305E40' and as a male parent the partially resistant line '67/3'. The fully resistant trait to Fom was introgressed from the allied species S. aethiopicum. In this work, the RIL population was assessed for the responses to Fom and by using a genomic mapping approach, two major QTLs on chromosomes CH02 and CH11 were identified, associated with the full and partial resistance trait to Fom, respectively. A targeted BSAseq procedure in which Illumina reads bulks of RILs grouped according to their resistance score was aligned to the appropriate reference genomes highlighted differentially enriched regions between resistant/susceptible progeny in the genomic regions underlying both QTLs. The characterization of such regions allowed us to identify the most reliable candidate genes for the two resistance traits. With the aim of revealing exclusive species-specific contigs and scaffolds inherited from the allied species and thus associated with the full resistance trait, a draft de-novo assembly of available Illumina sequences of the '305E40' parent was developed to better resolve the non-recombining genomic region on its CH02 carrying the introgressed Fom resistance locus from S. aethiopicum.

Ab initio GO-based mining for non-tandem-duplicated functional clusters in three model plant diploid genomes.

A functional Non-Tandem Duplicated Cluster (FNTDC) is a group of non-tandem-duplicated genes that are located closer than expected by mere chance and have a role in the same biological function. The identification of secondary-compounds-related FNTDC has gained increased interest in recent years, but little ab-initio attempts aiming to the identification of FNTDCs covering all biological functions, including primary metabolism compounds, have been carried out. We report an extensive FNTDC dataset accompanied by a detailed assessment on parameters used for genome scanning and their impact on FNTDC detection. We propose 70% identity and 70% alignment coverage as intermediate settings to exclude tandem duplicated genes and a dynamic scanning window of 24 genes. These settings were applied to rice, arabidopsis and grapevine genomes to call for FNTDCs. Besides the best-known secondary metabolism clusters, we identified many FNTDCs associated to primary metabolism ranging from macromolecules synthesis/editing, TOR signalling, ubiquitination, proton and electron transfer complexes. Using the intermediate FNTDC setting parameters (at P-value 1e-6), 130, 70 and 140 candidate FNTDCs were called in rice, arabidopsis and grapevine, respectively, and 20 to 30% of GO tags associated to called FNTDC were common among the 3 genomes. The datasets developed along with this work provide a rich framework for pinpointing candidate FNTDCs reflecting all GO-BP tags covering both primary and secondary metabolism with large macromolecular complexes/metabolons as the most represented FNTDCs. Noteworthy, several FNTDCs are tagged with GOs referring to organelle-targeted multi-enzyme complex, a finding that suggest the migration of endosymbiont gene chunks towards nuclei could be at the basis of these class of candidate FNTDCs. Most FNTDC appear to have evolved prior of genome duplication events. More than one-third of genes interspersed/adjacent to called FNTDCs lacked any functional annotation; however, their co-localization may provide hints towards a candidate biological role.

Segmental duplications are hot spots of copy number variants affecting barley gene content.

Copy number variants (CNVs) are pervasive in several animal and plant genomes and contribute to shaping genetic diversity. In barley, there is evidence that changes in gene copy number underlie important agronomic traits. The recently released reference sequence of barley represents a valuable genomic resource for unveiling the incidence of CNVs that affect gene content and for identifying sequence features associated with CNV formation. Using exome sequencing and read count data, we detected 16 605 deletions and duplications that affect barley gene content by surveying a diverse panel of 172 cultivars, 171 landraces, 22 wild relatives and other 32 uncategorized domesticated accessions. The quest for segmental duplications (SDs) in the reference sequence revealed many low-copy repeats, most of which overlap predicted coding sequences. Statistical analyses revealed that the incidence of CNVs increases significantly in SD-rich regions, indicating that these sequence elements act as hot spots for the formation of CNVs. The present study delivers a comprehensive genome-wide study of CNVs affecting barley gene content and implicates SDs in the molecular mechanisms that lead to the formation of this class of CNVs.

Transcriptomics, chromosome engineering and mapping identify a restorer-of-fertility region in the CMS wheat system msH1.

KEY MESSAGE: An original RNA-seq mapping strategy, validated with chromosome engineering and physical mapping, identifies candidate genes for fertility restoration in the 6HchS chromosome of Hordeum chilense in the wheat msH1 system. Cytoplasmic male sterility (CMS) is a valuable trait for hybrid seed production. The msH1 CMS system in common wheat results from the incompatibility between the nuclear genome of wheat and the cytoplasm of the wild barley Hordeum chilense. This work aims to identify H. chilense candidate genes for fertility restoration in the msH1 system with a multidisciplinary strategy based on chromosome engineering, differential expression analysis and genome mapping. Alloplasmic isogenic wheat lines differing for fertility, associated with the presence of an acrocentric chromosome Hchac resulting from the rearrangement of the short arms of H. chilense chromosomes 1Hch and 6Hch, were used for transcriptome sequencing. Two novel RNA-seq mapping approaches were designed and compared to identify differentially expressed genes of H. chilense associated with male fertility restoration. Minichromosomes (Hchmi), new smaller reorganizations of the Hchac also restoring fertility, were obtained and used to validate the candidate genes. This strategy was successful identifying a putative restorer-of-fertility region on 6HchS, with six candidate genes, including the ortholog of the barley restorer gene Rfm1. Additionally, transcriptomics gave preliminary insights on sterility and restoration networks showing the importance of energy supply, stress, protein metabolism and RNA processing.

Correction to: Identification and mapping of expressed genes associated with the 2DL QTL for fusarium head blight resistance in the wheat line Wuhan 1.

Following publication of the original article [1], we have been notified that some important information was omitted by the authors in the Copyright note. The Copyright note should read as below.

Grapevine comparative early transcriptomic profiling suggests that Flavescence dorée phytoplasma represses plant responses induced by vector feeding in susceptible varieties.

BACKGROUND: Flavescence dorée is the most serious grapevine yellows disease in Europe. It is caused by phytoplasmas which are transmitted from grapevine to grapevine by the leafhopper Scaphoideus titanus. Differences in susceptibility among grapevine varieties suggest the existence of specific genetic features associated with resistance to the phytoplasma and/or possibly with its vector. In this work, RNA-Seq was used to compare early transcriptional changes occurring during the three-trophic interaction between the phytoplasma, its vector and the grapevine, represented by two different cultivars, one very susceptible to the disease and the other scarcely susceptible. RESULTS: The comparative analysis of the constitutive transcriptomic profiles suggests the existence of passive defense strategies against the insect and/or the phytoplasma in the scarcely-susceptible cultivar. Moreover, the attack by the infective vector on the scarcely-susceptible variety prompted immediate and substantial transcriptomic changes that led to the rapid erection of further active defenses. On the other hand, in the most susceptible variety the response was delayed and mainly consisted of the induction of phytoalexin synthesis. Surprisingly, the jasmonic acid- and ethylene-mediated defense reactions, activated by the susceptible cultivar following FD-free insect feeding, were not detected in the presence of the phytoplasma-infected vector. CONCLUSIONS: The comparison of the transcriptomic response in two grapevine varieties with different levels of susceptibility to Flavescence dorèe highlighted both passive and active defense mechanisms against the vector and/or the pathogen in the scarcely-susceptible variety, as well as the capacity of the phytoplasmas to repress the defense reaction against the insect in the susceptible variety.

Genomic Regions From an Iranian Landrace Increase Kernel Size in Durum Wheat.

Kernel size and shape are important parameters determining the wheat profitability, being main determinants of yield and its technological quality. In this study, a segregating population of 118 recombinant inbred lines, derived from a cross between the Iranian durum landrace accession "Iran_249" and the Iranian durum cultivar "Zardak", was used to investigate durum wheat kernel morphology factors and their relationships with kernel weight, and to map the corresponding QTLs. A high density genetic map, based on wheat 90k iSelect Infinium SNP assay, comprising 6,195 markers, was developed and used to perform the QTL analysis for kernel length and width, traits related to kernel shape and weight, and heading date, using phenotypic data from three environments. Overall, a total of 31 different QTLs and 9 QTL interactions for kernel size, and 21 different QTLs and 5 QTL interactions for kernel shape were identified. The landrace Iran_249 contributed the allele with positive effect for most of the QTLs related to kernel length and kernel weight suggesting that the landrace might have considerable potential toward enhancing the existing gene pool for grain shape and size traits and for further yield improvement in wheat. The correlation among traits and co-localization of corresponding QTLs permitted to define 11 clusters suggesting causal relationships between simplest kernel size trait, like kernel length and width, and more complex secondary trait, like kernel shape and weight related traits. Lastly, the recent release of the T. durum reference genome sequence allowed to define the physical interval of our QTL/clusters and to hypothesize novel candidate genes inspecting the gene content of the genomic regions associated to target traits.

Identification and mapping of expressed genes associated with the 2DL QTL for fusarium head blight resistance in the wheat line Wuhan 1.

BACKGROUND: Fusarium head blight (FHB) is a problem of great concern in small grain cereals, especially wheat. A quantitative trait locus (QTL) for FHB resistance (FHB_SFI) located on the long arm of chromosome 2D in the spring wheat genotype Wuhan 1 is a resistance locus which has potential to improve the FHB resistance of bread wheat since it confers effective resistance to wheat breeding lines. Recently, differentially expressed genes (DEG) have been identified by comparing near isogenic lines (NIL) carrying the susceptible and resistant alleles for the 2DL QTL, using RNA-Seq. In the present study, we aimed to identify candidate genes located within the genetic interval for the 2DL QTL for FHB resistance, as assessed by single floret inoculation (FHB_SFI), and possibly contributing to it. RESULTS: Combining previous and additional bioinformatics analyses, 26 DEG that were located on chromosome arm 2DL were selected for further characterization of their expression profile by RT-qPCR. Seven of those DEG showed a consistent differential expression profile between either three pairs of near isogenic lines or other genotypes carrying the R and S alleles for the 2DL QTL for FHB resistance. UN25696, which was identified in previous expression work using microarray was also confirmed to have a differential expression pattern. Those eight candidate genes were further characterized in 85 lines of a double haploid mapping population derived from the cross Wuhan 1/Nyubai, the population where the 2DL QTL was originally identified. The expression QTL for gene Traes_2DL_179570792 overlapped completely with the mapping interval for the 2DL QTL for FHB_SFI while the expression QTL for UN25696 mapped near the QTL, but did not overlap with it. None of the other genes had a significant eQTL on chromosome 2DL. Higher expression of Traes_2DL_179570792 and UN25696 was associated with the resistant allele at that locus. CONCLUSIONS: Of the 26 DEG from the 2DL chromosome further characterized in this study, only two had an expression QTL located in or near the interval for the 2DL QTL. Traes_2DL_179570792 is the first expression marker identified as associated with the 2DL QTL.

Durum wheat genome highlights past domestication signatures and future improvement targets.

The domestication of wild emmer wheat led to the selection of modern durum wheat, grown mainly for pasta production. We describe the 10.45 gigabase (Gb) assembly of the genome of durum wheat cultivar Svevo. The assembly enabled genome-wide genetic diversity analyses revealing the changes imposed by thousands of years of empirical selection and breeding. Regions exhibiting strong signatures of genetic divergence associated with domestication and breeding were widespread in the genome with several major diversity losses in the pericentromeric regions. A locus on chromosome 5B carries a gene encoding a metal transporter (TdHMA3-B1) with a non-functional variant causing high accumulation of cadmium in grain. The high-cadmium allele, widespread among durum cultivars but undetected in wild emmer accessions, increased in frequency from domesticated emmer to modern durum wheat. The rapid cloning of TdHMA3-B1 rescues a wild beneficial allele and demonstrates the practical use of the Svevo genome for wheat improvement.

Transcriptomic and metabolomic analysis of ZmYUC1 mutant reveals the role of auxin during early endosperm formation in maize.

Kernel size in cereal is an important agronomic trait controlled by the interaction of genetic and environmental factors. The endosperm occupies most of the kernel area; for this reason, the endosperm cells dimension, number and metabolic content strongly influence kernel properties. This paper presents the transcriptomic and metabolomic analysis of the maize defective endosperm 18 (de18) mutant, where auxin accumulation in the endosperm is impaired. This mutation, involving the ZmYuc1 gene, leads to a reduced kernel size compared to the wild-type line B37. Our results mainly indicate that IAA concentration controls sugar and protein metabolism during kernel differentiation and it is necessary for BETL formation. Furthermore, a fine tuning of different auxin conjugates is reported as the main mechanism to counteract the auxin deficit. Some candidates as master regulators of endosperm transcriptional regulation mediated by auxin are found between MYB and MADS-box gene families. A link between auxin and storage protein accumulation is highlighted, suggesting that IAA directly or indirectly, through CK or ABA, regulates the transcription of zein coding genes. This study represents a move forward with respect to the current knowledge about the role of auxin during maize endosperm differentiation thus revealing the genes that are modulated by auxin and that control agronomic traits as kernel size and metabolic composition.

Omics approaches revealed how arbuscular mycorrhizal symbiosis enhances yield and resistance to leaf pathogen in wheat.

Besides improved mineral nutrition, plants colonised by arbuscular mycorrhizal (AM) fungi often display increased biomass and higher tolerance to biotic and abiotic stresses. Notwithstanding the global importance of wheat as an agricultural crop, its response to AM symbiosis has been poorly investigated. We focused on the role of an AM fungus on mineral nutrition of wheat, and on its potential protective effect against Xanthomonas translucens. To address these issues, phenotypical, molecular and metabolomic approaches were combined. Morphological observations highlighted that AM wheat plants displayed an increased biomass and grain yield, as well as a reduction in lesion area following pathogen infection. To elucidate the molecular mechanisms underlying the mycorrhizal phenotype, we investigated changes of transcripts and proteins in roots and leaves during the double (wheat-AM fungus) and tripartite (wheat-AM fungus-pathogen) interaction. Transcriptomic and proteomic profiling identified the main pathways involved in enhancing plant biomass, mineral nutrition and in promoting the bio-protective effect against the leaf pathogen. Mineral and amino acid contents in roots, leaves and seeds, and protein oxidation profiles in leaves, supported the omics data, providing new insight into the mechanisms exerted by AM symbiosis to confer stronger productivity and enhanced resistance to X. translucens in wheat.

Seed Dormancy Involves a Transcriptional Program That Supports Early Plastid Functionality during Imbibition.

Red rice fully dormant seeds do not germinate even under favorable germination conditions. In several species, including rice, seed dormancy can be removed by dry-afterripening (warm storage); thus, dormant and non-dormant seeds can be compared for the same genotype. A weedy (red) rice genotype with strong dormancy was used for mRNA expression profiling, by RNA-Seq, of dormant and non-dormant dehulled caryopses (here addressed as seeds) at two temperatures (30 °C and 10 °C) and two durations of incubation in water (8 h and 8 days). Aim of the study was to highlight the differences in the transcriptome of dormant and non-dormant imbibed seeds. Transcript data suggested important differences between these seeds (at least, as inferred by expression-based metabolism reconstruction): dry-afterripening seems to impose a respiratory impairment onto non-dormant seeds, thus glycolysis is deduced to be preferentially directed to alcoholic fermentation in non-dormant seeds but to alanine production in dormant ones; phosphoenolpyruvate carboxykinase, pyruvate phosphate dikinase and alanine aminotransferase pathways appear to have an important gluconeogenetic role associated with the restoration of plastid functions in the dormant seed following imbibition; correspondingly, co-expression analysis pointed out a commitment to guarantee plastid functionality in dormant seeds. At 8 h of imbibition, as inferred by gene expression, dormant seeds appear to preferentially use carbon and nitrogen resources for biosynthetic processes in the plastid, including starch and proanthocyanidins accumulation. Chromatin modification appears to be a possible mechanism involved in the transition from dormancy to germination. Non-dormant seeds show higher expression of genes related to cell wall modification, suggesting they prepare for acrospire/radicle elongation.

Comparative Transcriptome Profiles of Near-Isogenic Hexaploid Wheat Lines Differing for Effective Alleles at the 2DL FHB Resistance QTL.

Fusarium head blight (FHB), caused by the fungus Fusarium graminearum, represents one of the major wheat diseases worldwide, determining severe yield losses and reduction of grain quality due to the accumulation of mycotoxins. The molecular response associated with the wheat 2DL FHB resistance QTL was mined through a comprehensive transcriptomic analysis of the early response to F. graminearum infection, at 3 days post-inoculation, in spikelets and rachis. The analyses were conducted on two near isogenic lines (NILs) differing for the presence of the 2DL QTL (2-2618, resistant 2DL+ and 2-2890, susceptible null). The general response to fungal infection in terms of mRNAs accumulation trend was similar in both NILs, even though involving an higher number of DEGs in the susceptible NIL, and included down-regulation of the primary and energy metabolism, up-regulation of enzymes implicated in lignin and phenylpropanoid biosynthesis, activation of hormons biosynthesis and signal transduction pathways and genes involved in redox homeostasis and transcriptional regulation. The search for candidate genes with expression profiles associated with the 2DL QTL for FHB resistance led to the discovery of processes differentially modulated in the R and S NILs related to cell wall metabolism, sugar and JA signaling, signal reception and transduction, regulation of the redox status and transcription factors. Wheat FHB response-related miRNAs differentially regulated were also identified as putatively implicated in the superoxide dismutase activities and affecting genes regulating responses to biotic/abiotic stresses and auxin signaling. Altered gene expression was also observed for fungal non-codingRNAs. The putative targets of two of these were represented by the wheat gene WIR1A, involved in resistance response, and a gene encoding a jacalin-related lectin protein, which participate in biotic and abiotic stress response, supporting the presence of a cross-talk between the plant and the fungus.

Native soils with their microbiotas elicit a state of alert in tomato plants.

Several studies have investigated soil microbial biodiversity, but understanding of the mechanisms underlying plant responses to soil microbiota remains in its infancy. Here, we focused on tomato (Solanum lycopersicum), testing the hypothesis that plants grown on native soils display different responses to soil microbiotas. Using transcriptomics, proteomics, and biochemistry, we describe the responses of two tomato genotypes (susceptible or resistant to Fusarium oxysporum f. sp. lycopersici) grown on an artificial growth substrate and two native soils (conducive and suppressive to Fusarium). Native soils affected tomato responses by modulating pathways involved in responses to oxidative stress, phenol biosynthesis, lignin deposition, and innate immunity, particularly in the suppressive soil. In tomato plants grown on steam-disinfected soils, total phenols and lignin decreased significantly. The inoculation of a mycorrhizal fungus partly rescued this response locally and systemically. Plants inoculated with the fungal pathogen showed reduced disease symptoms in the resistant genotype in both soils, but the susceptible genotype was partially protected from the pathogen only when grown on the suppressive soil. The 'state of alert' detected in tomatoes reveals novel mechanisms operating in plants in native soils and the soil microbiota appears to be one of the drivers of these plant responses.

Corrigendum: Transcriptome Analysis of the Melon-Fusarium oxysporum f. sp. melonis Race 1.2 Pathosystem in Susceptible and Resistant Plants.

[This corrects the article on p. 362 in vol. 8, PMID: 28367157.].

Identification of bakanae disease resistance loci in japonica rice through genome wide association study.

BACKGROUND: Bakanae disease, caused by seed-borne Fusarium species, mainly F. fujikuroi, is a rice disease whose importance is considerably increasing in several rice growing countries, leading to incremental production losses. RESULTS: A germplasm collection of japonica rice was screened for F. fujikuroi resistance, allowing the identification of accessions with high-to-moderate levels of resistance to bakanae. A GWAS approach uncovered two genomic regions highly associated with the observed phenotypic variation for response to bakanae infection on the short arm of chromosome 1 (named as qBK1_628091) and on the long arm of chromosome 4 (named as qBK4_31750955). High levels of phenotypic resistance to bakanae were associated to the cumulated presence of the resistant alleles at the two resistance loci, suggesting that they can provide useful levels of disease protection in resistance breeding. A fine comparison with the genomic positions of qBK1_628091 and qBK4_31750955 with respect to the QTLs for bakanae resistance reported in the literature suggests that the resistant loci here described represent new genomic regions associated to F. fujikuroi resistance. A search for candidate genes with a putative role in bakanae resistance was conducted considering all the annotated genes and F. fujikuroi-related DEGs included in the two genomic regions highlighting several gene functions that could be involved in resistance, thus paving the way to the functional characterization of the resistance loci. CONCLUSIONS: New effective sources for bakanae resistance were identified on rice chromosomes 1 and 4 and tools for resistance breeding are provided.

Transcriptome Analysis of the Melon-Fusarium oxysporum f. sp. melonis Race 1.2 Pathosystem in Susceptible and Resistant Plants.

Fusarium oxysporum f. sp. melonis Snyd. & Hans race 1.2 (FOM1.2) is the most virulent and yield-limiting pathogen of melon (Cucumis melo L.) worldwide. Current information suggest that the resistance to race 1.2 is controlled by multiple recessive genes and strongly affected by the environment. RNA-Seq analysis was used to identify candidate resistance genes and to dissect the early molecular processes deployed during melon-FOM1.2 interaction in the resistant doubled haploid line NAD and in the susceptible genotype Charentais-T (CHT) at 24 and 48 h post-inoculation (hpi). The transcriptome analysis of the NAD-FOM1.2 interaction identified 2,461 and 821 differentially expressed genes (DEGs) at 24 hpi and at 48 hpi, respectively, while in susceptible combination CHT-FOM1.2, 882 and 2,237 DEGs were recovered at 24 hpi and at 48 hpi, respectively. The overall expression profile suggests a prompt activation of the defense responses in NAD due to its basal defense-related machinery that allows an early pathogen recognition. Gene Ontology (GO) enrichment analyses revealed a total of 57 GO terms shared by both genotypes and consistent with response to fungal infection. GO classes named "chitinase activity," "cellulase activity," "defense response, incompatible interaction," "auxin polar transport" emerged as major factors of resistance to FOM1.2. The data indicated that NAD reacts to FOM1.2 with a fine regulation of Ca2+-mediated signaling pathways, cell wall reorganization, and hormone crosstalk (jasmonate and ethylene, auxin and abscissic acid). Several unannotated transcripts were recovered providing a basis for a further exploration of the melon resistance genes. DEGs belonging to the FOM1.2 genome were also detected in planta as a resource for the identification of potential pathogenicity factors. This work provides a broader view of the dynamic changes of the melon transcriptome triggered by FOM1.2 and highlights that the resistance response of NAD is mainly signaled by jasmonic acid and ethylene pathways mediated by ABA and auxin. The role of candidate plant and fungal responsive genes involved in the resistance is discussed.

Comparative transcriptome profiling of resistant and susceptible rice genotypes in response to the seedborne pathogen Fusarium fujikuroi.

BACKGROUND: Fusarium fujikuroi is the causal agent of bakanae, the most significant seed-borne disease of rice. Molecular mechanisms regulating defence responses of rice towards this fungus are not yet fully known. To identify transcriptional mechanisms underpinning rice resistance, a RNA-seq comparative transcriptome profiling was conducted on infected seedlings of selected rice genotypes at one and three weeks post germination (wpg). RESULTS: Twelve rice genotypes were screened against bakanae disease leading to the identification of Selenio and Dorella as the most resistant and susceptible cultivars, respectively. Transcriptional changes were more appreciable at 3 wpg, suggesting that this infection stage is essential to study the resistance mechanisms: 3,119 DEGs were found in Selenio and 5,095 in Dorella. PR1, germin-like proteins, glycoside hydrolases, MAP kinases, and WRKY transcriptional factors were up-regulated in the resistant genotype upon infection with F. fujikuroi. Up-regulation of chitinases and down-regulation of MAP kinases and WRKY transcriptional factors were observed in the susceptible genotype. Gene ontology (GO) enrichment analyses detected in Selenio GO terms specific to response to F. fujikuroi: 'response to chitin', 'jasmonic acid biosynthetic process', and 'plant-type hypersensitive response', while Dorella activated different mechanisms, such as 'response to salicylic acid stimulus' and 'gibberellin metabolic process', which was in agreement with the production of gibberellin A3 in Dorella plants. CONCLUSIONS: RNA-seq profiling was performed for the first time to analyse response of rice to F. fujikuroi infection. Our findings allowed the identification of genes activated in one- and three- week-old rice seedlings of two genotypes infected with F. fujikuroi. Furthermore, we found the pathways involved in bakanae resistance, such as response to chitin, JA-dependent signalling and hypersensitive response. Collectively, this provides important information to elucidate the molecular and cellular processes occurring in rice during F. fujikuroi infection and to develop bakanae resistant rice germplasm.

Comparative Leaf and Root Transcriptomic Analysis of two Rice Japonica Cultivars Reveals Major Differences in the Root Early Response to Osmotic Stress.

BACKGROUND: Rice (Oryza sativa L.) is one of the most important crops cultivated in both tropical and temperate regions and is characterized by a low water-use efficiency and a high sensitivity to a water deficit, with yield reductions occurring at lower stress levels compared to most other crops. To identify genes and pathways involved in the tolerant response to dehydration, a powerful approach consists in the genome-wide analysis of stress-induced expression changes by comparing drought-tolerant and drought-sensitive genotypes. RESULTS: The physiological response to osmotic stress of 17 japonica rice genotypes was evaluated. A clear differentiation of the most tolerant and the most sensitive phenotypes was evident, especially after 24 and 48 h of treatment. Two genotypes, which were characterized by a contrasting response (tolerance/sensitivity) to the imposed stress, were selected. A parallel transcriptomic analysis was performed on roots and leaves of these two genotypes at 3 and 24 h of stress treatment. RNA-Sequencing data showed that the tolerant genotype Eurosis and the sensitive genotype Loto mainly differed in the early response to osmotic stress in roots. In particular, the tolerant genotype was characterized by a prompt regulation of genes related to chromatin, cytoskeleton and transmembrane transporters. Moreover, a differential expression of transcription factor-encoding genes, genes involved in hormone-mediate signalling and genes involved in the biosynthesis of lignin was observed between the two genotypes. CONCLUSIONS: Our results provide a transcriptomic characterization of the osmotic stress response in rice and identify several genes that may be important players in the tolerant response.